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Abstract. In equity and foreign exchange markets the risk- neutral dynamics of the underlying asset are commonly 
represented by stochastic volatility models with jumps. In this paper we consider a dense subclass of such models 
and develop analytically tractable formulae for the prices of a range of first-generation exotic derivatives. We provide 
closed form formulae for the Fourier transforms of vanilla and forward starting option prices as well as a formula for 
the slope of the implied volatility smile for large strikes. A simple explicit approximation formula for the variance swap 
price is given. The prices of volatility swaps and other volatility derivatives are given as a one-dimensional integral of 
an explicit function. Analytically tractable formulae for the Laplace transform (in maturity) of the double-no-touch 
options and the Fourier-Laplace transform (in strike and maturity) of the double knock-out call and put options are 
obtained. The proof of the latter formulae is based on extended matrix Wiener-Hopf factorisation results. We also 
provide convergence results. 
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1.1 Introduction 

A key step in the valuation and hedging of exotic derivatives in financial markets is to decompose these in 
terms of simpler securities, e.g. vanilla options, which trade in larger volumes, are generally very liquid and 
therefore have a well defined price. Such a decomposition is often achieved in two steps. First a model for 
the underlying asset under a risk neutral measure is calibrated to the implied volatility surface. In this step 
the current state of the market, as described by the prices of vanilla derivatives, is expressed in terms of the 
parameter values of the model. In other words the chosen model is used to impose a structure on the option 
prices. The second step consists of pricing the exotic derivative of interest in the calibrated model. 

It is well known that in equity and foreign exchange markets stochastic volatility models with jumps can 
be used to accomplish the first step described above (see e.g. [17], [24]). In the present paper we consider 
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forward starting vanilla options, volatility derivatives and barrier options, which are among the most widely 
traded exotic derivatives in the equity and foreign exchange markets. 

The desired properties of the model in each of the two steps described above place diametrically opposite 
restrictions on the choice of modelling framework. This is because in the calibration step one requires a 
flexible stochastic process that can describe well the current state of the vanilla market (i.e. can calibrate 
accurately to the observed implied volatility surface), while such flexibility can be a source of problems in 
the second step, where one needs to compute expectations of path-dependent functionals of the process. A 
more rigid modelling framework with, say, continuous trajectories and some distributional properties (e.g. 
independence of increments) could yield the structure of the process needed to establish efficient pricing 
algorithms for exotic derivatives. 

Wc investigate two families of stochastic volatility models with jumps: the time-changed exponential 
Levy models and the stochastic volatility models driven by Levy processes, where the volatility process is 
independent of the Levy driver. In these two families of models, the pricing of European derivatives is well 
understood and efficient calibration methods have been developed (see for example [12] and [24]), i.e. the first 
step of the two-stage procedure outlined above. However once the model is calibrated, the problem of pricing 
the first-generation exotic derivatives (e.g. barrier options) is quite involved. The law of the first-exit time 
from a bounded interval in stochastic volatility models with jumps, for instance, is not usually available in 
analytically tractable form. Because of the lack of structural properties that can be exploited to find the laws 
of the path-dependent functionals of interest, one would typically need to resort to Monte Carlo methods for 
the pricing of such derivatives in this setting. It is well-known that these methods are time-consuming and 
yield unstable results, especially when used to calculate the sensitivities of derivative securities. The method 
proposed in this paper to calculate the prices of such contracts consists of two steps: (i) a Markov chain 
approximation of the volatility process and (ii) an analytically tractable solution of the value function of 
the contract of interest in the approximating model. We provide proofs for the convergence of option prices 
under this approximation, and derive explicit expressions of Laplace and/or Fourier transforms of the value 
functions under the approximating model. 

The approximating class of stochastic volatility processes with jumps considered in this paper retains 
the structural properties required for the semi-analytic pricing (i.e. up to an integral transform) of forward 
starting options, volatility derivatives and barrier options. In the case of double-barrier option prices we 
will show that the process considered here admits explicit formulae for the Laplace/Fourier transforms in 
terms of the solutions of certain quadratic matrix equations. The main mathematical contribution of the 
present paper, which underpins the derivation of these closed form formulae, is the proof of the existence 
and uniqueness of the matrix Wiener- Hopf factorisation of a class of complex valued matrices related to the 
approximating model (see Theorem 4). These results extend those of [21] where the corresponding results 
for the real-valued case are established. In the context of noisy fluid flow models the matrix Wiener-Hopf 
factorisation for the case of a regime-switching Brownian motion is studied by [4, 28] . 
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It should be noted that the matrix Wiener-Hopf factorisation results developed in this paper can also 
be applied to the pricing of American call and put options in model (1-19) as follows. First we apply the 
main result in [21] to obtain the price of the perpetual American call or put and then, via the randomization 

algorithm introduced in [11], find the actual price of the option. 

Related Markov chain mixture models, which arc special cases of the model considered in this paper, have 
been studied before in the mathematical finance literature. In [19] and [22] explicit formulae were derived 
for the price of a perpetual American put option under a regime-switching Brownian motion model. The 
same process was used in [18] to model stochastic dividend rates where the problem of the pricing of barrier 
options on equity was considered. Finite maturity American put options were considered in [10] under a 
regime-switching Brownian motion model. More generally in [8], [7] numerical algorithms were developed 
in the case of regime-switching Levy processes. Furthermore extensive work has ham done on derivative 
pricing under stochastic volatility models with and without jumps (see the standard references [24], [17], [14] 
and [9]). 

The remainder of the paper is organized as follows. In Section 1.2 we state and prove the properties of 
continuous-time Markov chains and phase-type distributions that are needed to define the class of stochastic 
volatility models with jumps studied in this paper. In Section 1.3 we give a precise definition of this class 
of models and describe an explicit construction of an approximating sequence of models, based on Markov 
chains, which converges to the given stochastic volatility model with jumps. The models are set against the 
backdrop of a foreign exchange market which allows us to include naturally tlic stochastic foreign and domes- 
tic discount factors. In Section 1.4 we provide explicit formulae for the Fourier transforms in model (1.19) for 
vanilla and forward starting options. This section also gives an approximate explicit formula for the pricing 
of variance swaps, a one-dimensional integral representation of the price of a volatility swap and formulae 
for the asymptotic behaviour of the implied volatility smile for large strikes. Section 1.5 is devoted to the 
first-passage times of regime-switching processes. Section 1.6 discusses the pricing of double-no-touch and 
double-barrier knock-out options. It provides a formula for the single Laplace transform (in maturity) and 
the Laplace-Fourier transform (in strike and maturity) of the doublc-no-touch and the double-barrier knock- 
out options respectively in terms of the quantity that can be obtained from the complex matrix Wiencr-Hopf 
factorisation (see Theorem 6). Section 1.7 describes the fluid embedding of the model in (1.19), which plays 
a central role in the Wiener-Hopf factorisation. The key mathematical results of the paper, which allow 
us to price barrier options in the setting of stochastic volatility, are contained in Section 1.8 where matrix 
Wiener-Hopf factorisation is defined and the theorems asserting its uniqueness and existence are stated. The 
proofs of these (and other) results are contained in the appendix. 
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1.2 Mcirkov chains and phase-type distributions 



1.2.1 Finite state Markov chains 



We start by collecting some useful and well-known properties of finite state Markov chains that will be 
important in the sequel (see e.g. [16]). For completeness we will also present the proofs. Throughout the 
paper we will denote by 



M{i, j) = Mij = e'iMcj, m{j) = mj 



i,j = l,...,n, 



the ijih clement of an n x n matrix M and the jth clement of an n-dimcnsional vector m, where the vectors 
Ci, i = 1, . . . ,n, denote the standard basis of C" and where ' means transposition. Throughout the paper / 
will denote an identity matrix of appropriate size, and M+ = [0, oo) the non-negative real line. 

Lemma 1 Let Z he a Markov chain on a state space := {1, . . . , Nq], where Nq € N, and let B : E^ ^ C 
be any function. If Q denotes the generator of Z and Ab is a diagonal matrix of size Nq with diagonal 
elements equal to B{i), i = 1, . . . , Nq, then it holds that 



exp ^ B{Z,)dsj 
w/iereEi[-] =E[-|Zo = i], Pi[-] 



■Z. 



= e^Y>{t{Q + AB)){i,j) for any i,j€E'',t>0, (1.1) 
i], and I^.y is the indicator of the set {•}. 



Proof. Let (Pt)t>o be a family of A^o- dimensional square matrices with entries Pt{i,j), i,j = l,...,No, 
given by the left-hand side of (1.1). It is clear that Pq = I, where / is the A^o-dimensional identity matrix. 
The Markov property of the chain Z yields the Chapman-Kolmogorov equation Pt+s = PsPt — PtPs for all 
s,t>0. If we show that the family of matrices {Pt )t>o satisfies the system of ODEs with constant coefficients 



d P 

-^ = {Q + Ab)Pu 



Pa = L 



(1.2) 



then the lemma will follow, since equation (1.2) is well-known to have a unique solution given by the right- 
hand side of (1.1). The Chapman-Kolmogorov equation implies that Pt+h, — Pt = {Ph — I)Pt and it is therefore 
enough to show lim/i^o(-P/t — I)/h = Q + Ab- In other words we need to prove 



lim(P^(i,j)-/(i,i))//i 



Qihj) if«7^J, 
Qii,j)+B{j) ifi=j. 



(1.3) 



The random variables B{Zs) are bounded uniformly in s and hence the Taylor expansion of the exponential 
yields 



Ph{i,j)=Ei 



1 



B{Z,)ds 



+ o{h) for all i,j&{l,...No}. 



It is clear that 



1 Exotic derivatives under stochastic volatility models with jumps 



5 



lim I{z^=i\T 




B{Z,)ds = I{Zo=j}B{Zo) Pi-a.s. for all i,je{l,...No}, 



since the paths of Z are Pi-a.s constant for exponentially distributed amount of time. The dominated con- 



This concludes the proof. 

We now apply Lemma 1 to esatblish a simple but important property of the specturm of a discounted 
generator. 

Lemma 2 Let Q be a generator of a Markov chain with Nq G N states and let D be a complex diagonal 
matrix of dimension Nq . Then every eigenvalue X G C of the matrix Q — D (i.e. a solution of the equation 
{Q — D)x = Xx for some non-zero element x in C^°) satisfies the inequality 



where di = D{i, i), i = 1, . . . , Nq, are diagonal elements of D. In particular if min {^{di) : i = 1, . . . , Nq} > 
0, then the matrix Q — D is invertible. Furthermore, the real part of every eigenvalue of Q is non-positive. 

Proof. Let A be an eigenvalue of the matrix Q — D that corresponds to the eigenvector x G . Then x 
is also an eigenvector with eigenvalue exp(A) of the matrix exp((5 — D). Lemma 1 implies that if Z is the 
chain generated by Q then the following identity holds 



where (resp. di) denotes the i-th basis vector in C^" (resp. diagonal element of the matrix D). 

Assume now without loss of generality that the norm ||a;||oo '■= max{|a;j| : i = 1, . . . , Nq} of the vector x 
is one. Then identity (1.4) implies the estimate 



which proves the lemma. ^ 

Lemma 3 Let q £ C be such that 5R(g) > and M a matrix whose eigenvalues all have non-positive real 
part. Then the matrix ql — M is invertible and the following formula holds 



vergence theorem and the well-known fact Ej [-f{Zh=j}] = h{I{i,j) + Q{i,j)) + o{h) therefore imply (1.3). 



3?(A) < - min {^{di) : i = 1, . . . , iVo} 




(1.4) 



exp(SR(A)) = |exp(A)| = ||exp(Q-£))a;||oo <exp(-min{5R(di) : i = l,...,No}) 




(1.5) 



Proof. The following identity holds for any T G (0, oo) by the fundamental theorem of calculus 
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and, since the real part of the spectrum of the matrix M — ql is strictly negative by Lemma (2), in the limit 
as T ^ oo we obtain 

pOC 

exp((M - ql)t)dt = {ql - MyK 



I 

Jo 



□ 



1.2.2 (Double) phase-type distributions 

In this section we review basic properties of phase-type distributions, as these will play an important role in 
the sequel. Wc refer to Ncuts [27] and Asmussen [2] for further background on phase-type distributions. 

A distribution function F : IR.+ — > [0, 1] is called phase-type if it is a distribution of the absorption time 
of a continuous-time Markov chain on (m + 1) states, for some m e N, with one state absorbing and the 
remaining states transient. The distribution F is uniquely determined by the matrix A G M™^"*, which is 
the generator of the chain restricted to the transient states, and the initial distribution of the chain on the 
transient states a G K™ (i.e. the coordinates of a are non- negative and the inequalities < a'l < 1 hold, 
where 1 is the m-dimensional vector with each coordinate equal to one and ' denotes transposition). The 
notation X ^ PH(a,A) is commonly used for a random variable X with cumulative distribution function 
F. Note alse that the law of the original chain on the entire state space is given by 



the initial distribution and the generator matrix 





where denotes a row of m zeros. It is clear from this representation that the cumulative distribution 
function F and its density / are of the form 

F{t) = 1 - a'e*^l and f{t) = -a'e^^Al for any t G R+. (1.6) 

Note also that is an atom of the distribution if and only if a'l < 1 in which case the function / is a densitiy 
of a sub-probability measure on (0, oo). The n-th moment of the random variable X ~ PH{a, A) is given by 

E[X"] =n! a' (-^1)"" 1. 

It follows from the definition the phase-type distribution that the matrix A can be viewed as a generator 
of a killed continuous-time Markov chain on m states. Therefore we can express the matrix A as A = Q — D, 
where Q is the generator of a chain on m states and D \s & diagonal matrix with non-negative diagonal 

elements that are equal to the coordinates of the vector —Al. Lemma 2 therefore implies that the real part 
of each eigenvalue of A is non-positive. The next proposition gives a characterisation of the existence of 
exponential moments of a phase-type distribution in terms of the eigenvalues of the matrix A. 



Proposition 1 Let X ^ PH(a, A) be a phase-type random varaible as defined above and let Xq be the 
eigenvalue of the matrix A with the largest real part, i.e. SR(Ao) = max{5i(A) : A eigenvalue of A}. Then, for 
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any u G C, the exponential moment E[exp{uX)] exists and is finite if and only if^{u) < — 3?(Ao) in which 
case the following formula holds 

E[exp(uX)] = a' [A + uiy^Al + (1 - a'l), 
where I denotes an m- dimensional dentity matrix. 
Proof. It is clear that the identity 

POO 

E[exp{uX)] = F{X = 0) + / exp{tu)f{t)dt (1.7) 

Jo 

must hold for all w € M where / is the density of X on the interval (0, oo). Hence the question of existence of 
E[exp(?/X)] is equivalent to the question of convergence of the integral. Using Formula (1.6) for the density 
/, the fact cxp{t{A + ul)) = exp(tA) exp(to) for all u G C and the Jordan canonical decomposition of the 
matrix A we can conclude that 

E[| exp(uX)|] < 00 <^ -a' (^j^ exp((A + ^{u)I)t)dt^ Al< oo <^ SR(Ao + u) < 0, 

where Aq is as defined above. This proves the equivalence in the proposition. 

Note that the condition 3?(u) < — K(Ao) implies, by Lemma (2), that the matrix A + uI is invertible. For 
any T G R+ the fundamental theorem of calculus therefore yields the matrix identity 

/ exp{{A + ul)t)dt = {A + uI)-^[exp{{A + uI)T)-I]. (1.8) 
Jo 

Since all the eigenvalues of A + ul have a strictly negative real part, it follows from Jordan canonical 
decomposition of ^ + ul that limr^-oo exp{{A + uI)T) = 0. Therefore identities (1.7) and (1.8) conclude the 
proof of the proposition. ^ 

More generally, a double phase-type jump distribution DPH{p, I3~^ , B+ ,j3~,B~) is defined to have density 

/(x) :=p/+(.x)/(o,co)(a;) + (1 -p)/"(-a;)/(_oo,o)(a;) such that (1.9) 
pe[0,l], f^^PHi/3^,B^), /±(x) = -(/3±)V^*S±l and = 1, 



where the phase-type distributions PH{f3^ , B^) are as described above, 1 is a vector of the appropriate size 
with all coordinates equal to 1 and as usual Ia denotes the indicator of the set A. The condition l'/3^ = 1 
ensures that the distribution of jump sizes has no atom at zero. 

The class of double phase- type distributions is vast. Not only does it contain double exponential distri- 
butions 

/(x) := pa+e-^"^/(o,oo)(a;) + (1 -p)a"e'^""/(_oo,o)(a;) where > and p € [0, 1], (1.10) 

mixtures of double exponential distributions and Erlang distributions but this class is in fact dense in the 
sense of weak convergence in the space of all probability distributions on M. 
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Proposition 2 Let F be a probability distribution function on K. Then there exists a sequence (-Fn)neN of 
double-phase-type distributions Fn such that Fn=^ F as n — )• oo.^ 

This result directly follows from the three observations that (a) any probability distribution on the real line 
can be approximated in distribution arbitrarily closely by a random variable taking only finitely many values 
and (b) any constant random variable is the limit in distribution of Erlang or the negative of Erlang random 
variables, and (c) a mixture of Erlang distributions is a phase-type distribution. 

An important property of exponential distributions is the lack-of-memory property, which can be gener- 
alised to stopping times as follows: 

Lemma 4 Let {J^t)t>o « filtration and let p be any stopping timd^ with respect to this filtration. Let be 
an exponentially distributed random variable with parameter q > which is independent of the a -algebra 
generated by [Jt>oJ-'t- Then the equality 



E 



I{p<e^} exp(-A(eq - p)) 



^ e-iP holds for all A>0 



\ + q 



and hence the positive random variable Bg — p defined on the event {p < Bg} is, conditional on Tp, exponen- 
tially distributed with parameter q. 

Remarks, (i) This lemma can be viewed as a generalisation of the lack of memory property, 

P(eq > t + s\eg >s)= P(e, > t), 

of the exponential random variable when the constant time s is substituted by a stopping time p. Note 
also that it follows from the lemma that the conditional probability of the event {p < Sg} equals 

F{eg > p\Tp) = ex.p{-qp). 

(ii) Phase-type distributions enjoy a similar property that can be seen as a generalisation of the lack-of- 
memory of the exponential distribution. More specifically, let T follow a PH(a, B) distribution independent 
of the (T-algebra generated by Ut>oJ^t- Than for any stopping time p with respect to {J-'t)t>o, the random 
variable T — p defined on the event {p <T}, conditional on J^p, is PH(ap, B) distributed where 

Up = {a'eP^iy^a' exp{pB}, 

since the identity 



E 



J{^<T} exp(-A(r - p)) =a'exp{pB){B-XI) ^ Bl 



^ We write Fn =^ F for a sequence of distribution functions F„ and a distribution function F if Fn converges in 

distribution to F, that is, lim„_>oo Fn{x) = F{x) for all x where F is continuous. 
^ By definition the stopping time p takes values in [0, oo] and satisfies the condition {p < t} £ Ft for all t € [0, oo). 
The cr-algebra Fp consists of all events A such that Ar\{p <t} € Ft for all t £ [0, oo). 
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holds for all A > 0. This follows by the same argument as in the proof of Lemma 4. Furthermore we have 
the following expression for the conditional probability of the event {p <T}: 

P(T > plJ-p) = a'exp{pB)l. 



Proof. The following direct calculation based on Fubini's theorem, which is applicable since all the functions 
are non-negative, 



E [/A/{p<e,} exp(-A(eg - p))] = E 



POO 

J p 



E 



Ia 



Q 



X + q 



, where ^4 e J>, 



proves the identity in the lemma for all non-negative A. Since the Laplace transform uniquely determines 
the distribution of a random variable the lemma follows. ^ 



1.3 Stochastic volatility models with jumps 

We next describe in detail the two classes of stochastic volatility models with jumps that we will consider. 

Let V = {vt}t>o be a Markov process that takes positive values, modelling the underlying stochastic 
variance, and let X be a Levy process^ which drives the noise in the log-price process. The processes are 
taken to be mutually independent and are both defined on some probability space {f2,T,F). 

The law of X is determined by its characteristic exponent tp : M. M. which is according to the Levy- 
Khintchine formula given by 

E[e^"^*] = e*'^(") (1.11) 

with 

2 foo 

V'(u) = icu - —u^ + / [e^"" - 1 - iuxI^\,\<,}Hdx), (1.12) 

where cr^ > and c are constants and ly is the Levy measure that satisfies the integrability condition 
/ju(l A x'^)i'{dx) < oo. The triplet (c, cr^, i^) is also called the characteristic triplet of X. 

To guarantee that the option prices be finite we impose the usual restriction that X admits (positive) 
exponential moments; more precisely, we assume that for some p> 1 

/oo 
eP^^i^idx) < oo, (1.13) 

which implies that E[e^''^'] < oo for all t > 0. In this case the identity (1.11) remains valid for all u in the 
strip {u G C : ?s{u) G {—p,0]} in the complex plan, where the function ijj is analytically extended to this 
strip. 



® A Levy process X = {.X't}t>o is a stochastic process that has independent and stationary increments, and has 
right-continuous paths with left-limits with Xo = 0. 
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In this setting a candidate stochastic volatiUty process with jumps S = {St}te[o,T]j where T > denotes 
a maturity or time-horizon, is given by 

:= So exp ^(r - d)t + \/^dX„ - ^ ^(-iv^)ds^ , Sq = s > 0, (1.14) 

where r and d are the instantaneous interest rate and dividend yield respectively. Here we assume that the 
variance process v satisfies the following integrability condition: 

i-T 

\tp{—i^yv^)\ds < 00 a.s. (1-15) 



/ 

Jo 
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It is easy to see by conditioning on the filtration generated by the variance process v that the integrability 
condition in (1.15) implies the martingale property of the discounted process {c~^'^~'^^*St}te[a.T]- 

Note that if wc take for example X to be a Brownian motion with drift and v an independent square- 
root process, the process S reduces to a Hcston model with zero correlation between the driving Brownian 
factors (see e.g. [17]). The class of models described by (1.14) is quite flexible, and contains for example the 
stochastic volatility models with jumps described in Lipton [24], as long as there is no correlation between 
the driving Brownian motions. 

A related class of models that has been proposed in the literature is the one where the effect of stochasticity 
of volatility is achieved by randomly changing the time-scale (see e.g. Carr et al. [12]); in the setting above 
the price process {St}te[o,T] is defined by 

St := exp ((r - d)t + Xy, - ip{-±)Vt) , 5*0 = s > 0, where Vt := [ Vudu, (1.16) 
and we assume that v satisfies the integrability condition 

Vt < oo a.s. (1.17) 

Also in this case the discounted process {e~^'^~'^^*St}te[o^T] is a martingale. 

It is clear from the definitions that in the case where X is a Brownian motion with drift, the classes of 
models in (1.14) and (1.16) coincide, due to the scaling property of Brownian motion. Whereas the effect of 
the variance process v on the Brownian motion with drift is the same in both classes of models, the effect 
of the process v on the behaviour of jumps is different. In (1.16) the Markov process v modulates only the 
intensity of the jumps of X while in model (1.14) the volatility scales the distribution of size of the jumps 
but does not affect the intensity. 

In the next section we will describe a modelling framework in which any model in the classes given 
by (1.14) and (1.16) can be approximated. The approximation in Section 1.3.2 retains the structural prop- 
erties required for the semi-analytic pricing (i.e. up to an integral transform) of barrier options, forward 
starting options and volatility derivatives. 
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1.3.1 A class of regime-switching models 



Let the set := {1, . . . , Nq} be the state-space of a continuous-time Markov chain Z = (.^f)t>o and let the 
the process W = {Wt)t>o denote a standard Brownian motion which is independent of the chain Z. For each 
i G let the process J' := {Jl)t>o be a compound Poisson process with intensity Aj > and jump-sizes 
distributed according to a double-phase- type distribution DPH{pi , , , /3r , ) . In particular the jump- 
size distributions have no atom at zero, i.e. (13^)' 1 = 1 for all i G E^ such that A,p, > and analogously 
for Assume further that the processes J' are mutually independent as well as independent from the 
Brownian motion W and the chain Z. 

In this setting consider the following model for the underlying price process S = (S'j)t>o, the (domestic) 
money market account = (i3^)t>o and the cumulative dividend yield B^ = {B[)i>Q: 

Bf := exp (^J^ RD{Zs)ds^ , Sf := exp (^J^ RF{Zs)ds^ , St := exp(Xt), (1.18) 



where 



Xt:=x + [ fiiZs)ds + [ a{Z,)dW, + I hz.=i}^Ji- (1-19) 

In the case of the Foreign Exchange market the process B^ can be interpreted as a foreign money market 
account. The point x G R is the starting value of the process X and Rd,Rf, ^J^,cr : -E" — )■ R are given real- 
valued functions on E^ such that Rd,Rf are non-negative and ct is strictly postitive. To price derivatives 
in our model, we need to understand the law of the Markov process {X,Z), which is determined by the 
characteristic matrix exponent K, defined as follows. 

Definition. The characteristic matrix exponent K -.'^^ (^NqxNo ^j^", Z) is given by 

K{u) = Q + A{u), 

where Q denotes the generator of the chain Z and, for u S K, A{u) is a diagonal matrix of size No x Nq, 
where the i-th diagonal element equals the characteristic exponent of the process X in regime i, given by 

i^iiu) := iwrt - a^V2 + Ai [pi(/3,^)'(B+ + iuiy^B+l + (1 - p,)(/3r)'(Br - iuiy'Bri - l] . (1.20) 

where / and 1 are an identity matrix and a vector with all coordinates equal to one of the appropriate 
dimensions. 

Remarks, (i) Note that the fimctions Vi defined in (1.20) can be analytically extended to the strip in the 
complex plane ?s{u) G (— Q;^,a^) where 

af := min{— 5ft(A) : A eigenvalue of B^} for any state i G E'". (1-21) 



(ii) In the special case that the jumps follow a double exponential distribution the diagonal elements of 
the matrix A{u) take the simpler form 
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= ±uni - a>V2 + ( . - + Ai(l - Pi) ^ 



where and Pi are the parameters of the double exponential distribution. 

(iii) Throughout the paper we will use Ea:^j["] to denote the conditional expectation E[-|Xo = x,Zo = i] 
and on occasion to represent Eo,i[']- 

We now define two matrices the will play an important role in the sequel. 
Definition. The discount rate matrix Ad is the diagonal matrix with elements A£){i,i) '■= Roii) where 
i € E'^. The dividend yield matrix Ap is the diagonal matrix given by AF{i,i) := -RfW for i G 

Theorem 1 The discounted characteristic function of the Markov process {X, Z) is given by the formula 

exp(iuXt) 



"-^v J =exp{±ux) ■exp{t{K{u) - AD)){i,j) (1.22) 
L -D* J 

for all ueR. 

Remarks, (i) The left-hand side is finite for all w € C in the strip G (— a+,a~) where 

at = min{a^ : XkPk >0,k& E°} and a' = mm{a^ : Afe(l - pk) >0,ke E°}, (1.23) 

the quantities af are defined in (1.21) and the minimum over the empty set is taken to be +oo. It follows, 
by analytical continuation, that the identity (1.22) remains valid for all u in this strip. Furthermore, the 
slope of the implied volatility smile in model (1-19) is determined by a+ and a~ (see Subsection 1.4.2). 

(ii) The Markov property and Theorem 1 imply that the process {SfBf /B^}t>o is a martingale if the 
following two conditions hold: 

1 < a^, for all k G E° such that Xpk > 0, (1.24) 
A{-i) = Ad- Ap. (1.25) 

Condition (1.24) ensures that Ei^2;[5T] is finite for all T > and hence by Theorem 1 takes the form 
Ei,a;['S'T] — [exp {TK{—i)) 1] {i). The equality in (1.25) guarantees that S has instantaneous drift given by 
the rates Ad — Ap. Any model from the class (1.18)-(1.19) that satisfies conditions (1.24) and (1.25) can 
be taken as a specification of the price process of the risky-asset under a pricing measure. Prom now on we 
assume that model (1.19) is specified under the pricing measure given by condition (1.24)-(1.25). 

(iii) For later reference we record that, under a pricing measure, the price at time s of a zero coupon 
bond maturing at time t > s is given by 



Ei 



5f 



77(^,2) 



• (exp((i - s)iQ - Ad))1) {Zs), (1.26) 



where J^^^''^^ = a{{Xu, Zu)}u<s denotes the standard filtration generated by {X, Z). In particular, at time 
the price is given by 
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(Bf ) '1 = (exp(t(Q - An))l) {i). 



(iv) The infinitesimal generator £ of the Markov process {X, Z) acts on sufficiently smooth functions^ 
/ : R X R as 



Cf{,x, i) = (x, i) + i) + 



/ f{x-\- z,i)gi{z)dz- f{x,i) 



+ E (1-27) 

where gi is the density of the double phase- type distribution DPH{pi, , , ,B^), qij is the ij th 

element of Q and ' denotes differentiation with respect to x. 

(v) For a specific regime-switching model (namely the case where the Markov chain Z has two states 
only) the calibration is studied in [26]. 

Proof. It is clear from the definition of {X, Z) that it is a Markov process. Let J'f := (j{Zs : s G [0,t]) be 
the fT-algebra generated by the chain Z up to time t. Since the compound Poisson processes and Brownian 
motion in model (f .19) are mutually independent as well as independent of the cr-algebra , it is easy to 
see that by conditioning on J^f for any i e we obtain 

E^,i [exp{uXt)\J^f] = exp (ux + uj^ n{Z,)ds + y ^ (riZ^fds + v{Zs, u)d^ , (1.28) 
L'{i,u) := Xi [E[exp(uJj)] - 1] 

= A, [MPtYiBt + uir'B+l + (1 - Vim)'{Br - uiy^Bri - l] , 

where the random variable Ji denotes the size of jumps of the compound Poisson process J*. The last equality 
in this calculation is a consequence of the choice (1.9) of the distribution of jump sizes and Proposition 1. 
Therefore the complex number u must be contained in all intervals {—a'^,a'^), k e E^, where are defined 
in Theorem 1. The identity in (1.28) holds more generally for any jump-distribution that admits a moment 
generating function. The well-known identity from the theory of Markov chains given in Lemma 1 can now 
be applied to obtain the expectations of the expressions on both sides of (1.28). This concludes the proof of 
Theorem 1. 

□ 

Prom Theorem 1 one may obtain an explicit expression for the marginal distributions of {X, Z) by 
inverting the Fourier transform (1.22): 

Proposition 3 For any T > 0, the joint distribution q^^iy^j) = ■^fx,i[XT < y^Zx = j] of {Xt,Zt) is 
given by 

qT'\y,j) = ^ [ e'^^^-y^e^p{K{^)T){i,j)d^, yGR,i,jGE°, (1.29) 



® For example, functions / with f{-,i) € C^{^) for i € E°, where C^{R) are the twice continuously differentiable 
functions with compact support. 
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In particular, Xt is a continuous random variable with probability density function q^'^{y) = ''"''''^^^^'^^^ 
given by 

iViy) = ^ [ [exp {K (0 T) 1] (z) d^, y e M, z e (1.30) 

Proof. It is well-known that a probability law on the real line M has a density with respect to the Lebesgue 
measure if its characteristic function is in i^(M). The characteristic function of Xt by Theorem 1 equals 

E,,,- [exp(ieXT)] = e^«-'^-[cxp(i^(OT)l](i). (1.31) 

We now show that this characteristic function is asymptotically equal to cxp(— c^^), as |^| — >■ oo, for some 
positive constant c. Note first that the volatility vector a in model (1.19) has non-zero coordinates by 
assumption and the spectra of matrices , i = 1, . . . ,No, do not contain any points of the form i^, for 
^ e M, by Lemma 2. Therefore the functions ^ i-)- i^i{£,), i = 1, . . . ,iVo, defined in (1.20) are asymptotically 
equal to downward facing parabolas. A further application of Lemma 2 implies that the characteristic function 
has the desired asymptotic behaviour. This further implies that the density g^'* exists and is given by the 
inversion formula (1.30). 

Theorem 1 and a similar argument to the one outlined in the previous paragraph imply that the function 

^ ^ e-i- eMTK{0){i,j) = [exp(i^XT)/{z.=,}] 

is in L^(]R) and that the two equalities hold. Therefore the Fourier inversion formula is valid and the identity 
in (1.29) follows. ^ 



1.3.2 Two step approximation procedure 

The construction of a regime-switching Levy process with jump sizes distributed according to a double 
phase- type distribution that approximates a given stochastic volatility process with jumps from either of the 
two-classes (1.14) and (1.16) takes place in two steps: 

(i) Approximation of the variance process w by a finite-state continuous-time Markov chain and 

(ii) Approximation of the Levy process X by a Levy process with double-phase-type jumps. 

By approximating the variance process by a finite state Markov chain the resulting approximating process is 
a regime-switching Levy process. The approximation of the jump part of X by a compound Poisson process 
with double phase- type jumps will enable us to employ matrix Wiener-Hopf factorisation results, needed to 
obtain tractable formulae for the prices of barrier-type options. The two steps will be described in detail in 
the present section. 
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Markov chain approximation of the variance process 



The first step of the approximation procedure that was outlined above is to approximate the variance process 
w by a finite-state continuous-time Markov chain on some grid contained in the positive real line. We will 
restrict ourselves to the case that the variance process w is a Feller process on the state space K+ = [0, oo). 
This assumption implies that w is a Markov process satisfying some regularity properties. 

The Feller property is phrased in terms of the semi-group {Pt)t>o of v acting on Co(IR+), the set of 
continuous functions on 1R+ that tend to zero at infinity. Recall that, for any Borel function / on R_|_ and 
t>0, the map Ptf : M+ -J> R is given by 

PJix)=EM{vt)]. 

Assumption 1 The Markov process v = {vt}t>o is a Feller process; that is, for any f G Co(M+), the family 
{Ptf)t>o satisfies the following two properties: 

(i) Ptf e Co(K+) for any t > 0; 

(ii) limtio Ptf{v) = f{v) for any v gR+. 

An approximating Markov chain Z with generator Q on a state-space E° = {xi,. . . , xnq} can be con- 
structed by choosing to be some appropriate (non-uniform) grid in M+, and specifying the generator Q 
such that an appropriate set of instantaneous (local) moments of the chain Z and the target process v are 
matched. Sec [25] for details on this procedure. 

Denote hy Q : V ^ Co(M+) the infinitesimal generator of v defined on its domain V, and let Z'^"' be a 
sequence of Markov chains with generators Q^"^ and state-spaces 

^o(n) ^ ^ _ _ ^ x'-^l^}, and denote by 

Q^"'^fn the vector with coordinates 



Q^"^fn{Xi)= Q^"\xi,Xj)f{Xj), Xi€E^^^\ 

The sequence Z^") weakly approximates the variance process v if the range of the state-spaces E^^"'^ grows 
sufficiently fast as n tends to infinity, and if, for all regular functions /, Q^"'^fn converges uniformly to Qf, 
that is e„(/) — )• where 

e„(/):= max Q^""^ fn{x) - G f (x) 

a;e (£;"("')" 

and (£;"(") )° is equal to E^^"'^ without the smallest and the largest elements. The precise statement reads as 
follows: 

Theorem 2 Assume that the following two conditions are satisfied for any function in a core' of C: 

^nif) — > as n ^ oo, (1-32) 
either (i) lim gfiy) = or (ii) lim P^^ [r,^"Lu„ > = 1, (1.33) 



A core C of the operator £ is a subspace of the domain of £ that is (i) dense in Co(K+) and (ii) there exists A > 
such that the set {(A — £)/ : f € C} is dense in Co(]R+). 
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where for any set G C K+ we define Tq^ = inf{t > : Xj"^ ^ G}. Then it holds that, as n ^ oo, t;'"^ ^ v.^ 
Remark, (i) Note that the convergence in law implies in particular that 



(4"))] ^E,\g{vT)] 



for any bounded continuous function g. 

(ii) A proof of this statement can be found in [25]. 

Approximation of a Levy process 

The second stage of the aforementioned approximation procedure amounts to an approximation of a Levy 
process by a compound Poisson process with double phase-type jumps. 

Proposition 4 For any Levy process X there exists a sequence of Levy processes with double 

phase-type jumps such that X^"^ ^ X as n —> oo. 

Remarks, (i) A proof of this result can be found in e.g. Jacod and Shiryacv [20, Section VII. 3]. It is based 
on the fact that a sequence (X("^)„gN of Levy processes weakly converges to a given Levy process X if and 
only if converges in distribution to Xi (see e.g. [20, Corollary VII. 3. 6] for a proof). 

(ii) Furthermore [20, Corollary VII. 3. 6] implies that a sufficient condition to guarantee that x["'^ converges 
in distribution to Xi is that the characteristic triplets (c„,cr^,f„) of X'"^ converges to the triplet (c, cr^,i/) 
of X as follows as n — >■ oo : for some a > that is a continuity point of ^{dx) and ^{—dx) it holds that 

c„ -)> c, o-^ + / x^fc'„(da;) -J> cr^ + / x'^i'{dx) (1-34) 

J {—a, a) J (—a, a) 

/ (a;^ Aa)|F„(a;) -F(a;)|da; + / (|a;|^ A a)|i/„(a;) - z/(a;)|da; (1.35) 

J(0,oo) i(-(x>,0) 

where, for any measure m on M, m and m are the left and right tails, rri{x) = m([a;, oo)), m{x) = m((— oo, x\). 

Suppose now that the Levy process X is a model for the log of a stock price. Then X has a triplet (c, , v) 
satisfying the exponential moment condition in (1.13). An example of a sequence of Levy processes with DPH 
distributed jumps that weakly converges to X is then given as follows. Let A„ = y[{—l/n, 1/nY) and F„ 
be a double phase- type distribution that approximates in distribution the probability measure Fn{dx) = 
I{\x\>i/n}^{'^x)/\n, and define the measure Vn by yn{dx) = A„F„(da;). Here the F„ and are to be chosen 
such that (1.35) and the second requirement in (1.34) hold true. See e.g. [3] for a fitting procedure based 
on the EM algorithm. Then the sequence of Levy processes A"'"' with triplets (c, cr^,A„F„) satisfies the 
conditions (1.35) and thus approximates X in law as n tends to infinity. 



® By ^ v we denote the convergence in law of u'"^ to v in the Skorokhod topology, i.e. convergence of the 
distributions of u*-"^ to those of v in the set of probabiUty measures on the Skorokhod space endowed with the 
Skorokhod topology. 
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1.3.3 Convergence of the approximation procedure 



Combining the two steps in the approximation we can now identify candidate sequences of regime-switching 
processes tliat converge to either of the stochastic volatihty processes with jumps (1-14) and (1-16) and 
establish the convergence. 

Let (X("))„eN be a sequence of Levy processes with DPH jumps, and let (Z("))„eN be a sequence of 
Markov chains that is independent of X and Let ipn denote the characteristic exponent of X^"), 

and (c„,(T^,t'„) the characteristic triplet. Consider then the sequences of stochastic processes 
and (5(*--"))„eN with = {St'~''^}te[o,T] and = {st~''^}t^io,T] given by 

giint-n) _ ^ y^dXi") - £ (^i^^) ds^ , 

5f := So exp ((r - d)t + (v/")) - ^„(-i)V^/")) , 

pt 

where v/"^ = / ^("Ms. 

Jo 

The processes S'C*"*-") and S^*""") are in law equal to exponential Levy processes: 

Proposition 5 For n G N, logS'^'"*"") and logS^*^"") are in law equal to regime- switching Levy processes 
of the form (1.19). 

Proof. Note first that since X^"^ = (Xj"^)t>o is a Levy process with DPH jumps it is of the form 

^ = ii^^h + a^^'^Wt + Ji"^ with J,("' = Ui, (1.36) 

where ij,^"'\a^"'^ are constants, M'"^ = (M("^)(>o Poisson processes with jump-rates A'"^ and Ui are i.i.d. 
random variables following a DPH distribution. Then it is clear that X^'"*"") = log(S'^'"*~"V'5'o) is in law 
equal to the process x^*"*"") = (X(*"*~"^)t>o given by 

Aft") t 

where J("'^),j = 1, . . . , A/'^"), are independent compound Poisson processes, that are in law equal to the 
processes = {Jt^'^^)t>o with J^""'^^ = J2f=i \J x^^Ui, respectively. Since, for any constant c 7^ 0, cUi 

follows a DPH distribution, is a regime-switching Levy process of the form (1.19). 

As a consequence of the scaling property of Brownian motion it follows that = log(5(*'^~"^/S'o) 

is in law equal to the process = (X(*'^~"^)t>o given by 
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. t 

(tc—n) 



where J("'^),j = 1, . . . , A/'f"), are independent compound Poisson processes, that are in law equal to the 
processes = (Jt"'''^)t>o with J^"'-'^ = f/j, respectively, where M^"'^^ is a Poisson process with 

jump-rate A • . Here we used that, conditional on Z^"\ the process X^*"^""^ has independent increments, 
so that the law of X^**^""^ conditional on Z"-"' is determined by the conditional characteristic functions of 
^{tc ^ > Q_ ^ straightforward calculation verifies that the conditional characteristic functions of x'f" 
and xi*'^""^ are equal for t> 0. „ 

To establish the convergence in law of (5('"*-"))„eN and (S'(*'=-"))„eN we will first show convergence of 
the finite-dimensional distributions.^ 

Proposition 6 Assume that X and V as n tends to infinity. Then the following holds true: 

(a) (Z("),5('"*-"))„ {v,S) where S is the model given in (1.14). 

(b) 6'(*'^~"^)„ {V,S) where S is the time-change model given in (1.16). 

Remark. The convergence of European put option prices (and hence, by put-call parity, also of European 
call option prices) under the approximating models to those under the limiting models is a direct consequence 
of the convergence in finite dimensional distributions. To establish the convergence of path-dependent option 
prices such as barrier option prices it is required to prove that the approximating models converge in law. 

Proof, (b) To prove the convergence of the finite dimensional distributions it sufltices, in view of the Markov 
property, to show that, for each fixed t e [0, T], the characteristic functions Xn of (v/"\ X(")(V^*'"'')) converge 
point-wise to the characteristic function x of {yt,X{Vt)) as n tends to infinity. By conditioning and using 
the independence of V'"' from X^") and of V from X we find that 

Xn{u,v) = E^^i exp|(iu-h Vn(^^))v/"'^| , x{u,v) =E^^i[exp{{±u + ip{v))Vt)]. 

Since v^"'^ converges in law to v in the Skorokhod topology, the Skorokhod representation theorem implies 
that on some probability space Z^"'^ — )• v, almost surely, with the convergence with respect to the Skorokhod 

metric. Since, for any t e [0,T], the map it : -DR[0,r] — >■ R given by li : x i— )• J*x{s)ds, is continuous in the 

(n) 

Skorokhod topology, we deduce that — >■ Vt almost surely. In particular, Xn converges point-wise to x- 
The proof of (a) is similar and omitted. ^ 

The next result concerns the convergence in law of the sequences (^S^^^^ ^^)?ign 

and (S'(*^-"))„eN: 



A sequence of processes (yn)neN converges in finite dimensional distribution to tiie process Y = {^^t}te[o,T], if, for 
any partition ti < . . . < tm of [0, T], ¥{y}"^ <xi,..., y}^^ < Xm) P(Fti <xi,...,Yt^ < Xm). We will denote 
this convergence by Yn Y. 
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Theorem 3 The following statements hold true: 

(a) Assume that {v,S) is a Feller process, where S is the model given in (1.14), that Z^") satisfies the 
conditions in Theorem 2, and that the characteristics of X^'^^ satisfy conditions (1-34) and (1.35). Then, as 

oo, 

(b) Assume that 4> X and Z^"^ 4> w as n ^ oo. Then it holds that 

S(tc-n) 4 

where S is the time-change model given in (1.16). 

Remark. The convergence in law stated above carries over to the convergence of barrier option prices under 
the respective models, if the boundaries are continuity points of the limiting model. For instance, itwc denote 
by TA = inf{t > : St <^ A} the first time that S leaves the set A := [l,u], and P{St e {i,u}) = 0, then, 
for any bounded continuous pay-off functions g,h : R+ — >■ R, we have that as n — oo the double knock-out 
option and rebate option prices under the approximating models converge to those under the limiting model: 



>T}\ [9{ST)I{rA>T}] 



{rr<T} 



[e— /i(5,J7{,,<T}] , 



where 5^") denotes S^*"* or S^*" . A proof of this result was given in [25] . 

Proof. In view of Proposition 6, it sufBces^^ to verify that the sequences (S'(*"*~"))„£n and (<S'^*'^~"^)neN are 
relatively compact in Dk[0,T]. 

(a) Wc will establish relative compactness of the sequence (X(*"*~"^)„gN = (log S'('"*^"^)„eN. Let X' = 
log S. It is straightforward to check that the set of functions / : M+ x R+ — > M of the form f{x, v) = g{x)h{v) 
with h in the core of Q, the infinitesimal generator of v, and with g G C^(IR)^^ is dense in Co(M+ x R) and 
is contained in the domain 'D{jC') of the infinitesimal generator C of {v,X'). Furthermore, £' acts on such 
/ as 

^'fix, v) = ^vhiv)a^g"ix) + [{r - d) + - V;(-iv^)] h{v)g'{x) 



+ 



K"") I [9{x + z^/v)-g{x)-z^/vg\x)I{\^\<l}\v{dz)-\-g{x)gh{v), xeR,w>0, 



since by construction the stochastic integral |jJy^dXs| jumps if and only if the Levy process X 
jumps and, if the jump occurs at time t, the quotient of the jump sizes equals y^. On the other hand, the 
infinitesimal generator of the regime-switching processes (Z^"^ X^'"*""^) acts on f{x,v) = g{x)h{v) as 



See e.g. Theorem 3.7.8 in Ethier and Kurtz [15] for a proof of this well known fact. 

C^(R+) denotes the set of infinitely differentiable functions with compact support contained in K+. 
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v) = ^valh(v)g'\x) + [{r - d) + c„V^ - ^„(-iV^)] h(v)g'{x) 

+ h{v) I [g{x + zV^)-g{x)-zV^g'{x)I{\,\<iy]un{dz)+g{x)Q^"'^h{v), x€R,v€E°^''\ 
Js. 

C^"'^f converges to Cf uniformly as n ^ oo: 

e;(/):= sup \£'f{x,v)-£^"'>f{x,v)\^0. (1.37) 

To see why this is true note that the triangle inequality and integration by parts imply that 

e'M) <\\Gh~ Q(")/i||„||,g||co + C\\v{a) - I7„(a)| + C2 max{\al - a% \cn - c|} 
+ C3|^"|F(^)-F„(^)|d^ + j \v{z)-v^{z)\dz^ 

+ cA[ z^\V{z)-Uniz)\dz+ [ z^\ii{z) - i^Jz)\dz\ , (1.38) 

lJ{0,a) J{-afi) J 

where a is a continuity point of the measures u{dx) and i'{—dx) and Ci, . . . , C4, are certain finite constants 
independent of n, and we denoted ||/||„ = sup3.£^o(n) \f{x)\ and ||/|| 00 — ^uPxGM view of the 

conditions (1.34) and (1.35) e'j(/) tends to zero as n tends to infinity. Corollary 4.8.6 in Ethier and Kurtz [15] 
implies then that (Z("), X(™*-"))„eN and hence (Z("), 5(™*""))„eN is relatively compact in D^2[0,T]. 
(b) Denote by X^"' — {Xf"^}t>o and X = {X^jt^Q the Levy processes given by 

^in) ^ -^(n) _ ^nd Xt = Xt - i^{-±)t. 

We will verify^^ the relative compactness of the sequence {Y'^"^)nQff with Y^"'^ = X^"'\Vf "'^). In view of 
the Skorokhod embedding theorem and the convergence in law of (^("^X^")) to {v,X), we may and shall 
assume that (Z("\ X("^)„gN and {v,X) are defined on the same probability space (i?', J^',P'), and that, 
P'-almost surely, (Z("\X(")) {v,X) with respect to the Skorokhod metric. Fix an a; e /?' for which this 

convergence holds true. 

Observe that, for any U > 0, z = sup„ Z^^\uj) is finite, as vu^cu) is finite and \Z^\uj) — vu{u}\ — > as n 
tends to infinity. For x G -De[0, U], 6 > 0, U > 0, denote by 'w'{x, 6, U) the modulus of continuity 

w'{x,S,U) = inf max sup \x{s) — x{t)\, 
» s,te[ti-i,ti) 

where {ti} ranges over all partitions = to < ti < . . . < tn-i <U <tn with mini<j<„ \ti — > 5. Note 
that w'{x, 5, U) is non-decreasing in 5 and U. Therefore it is straightforward to check that 

w'^{U) :=w;'(y(")(w),(5,C/) < w' {X^''\w),5z,Uz). (1.39) 

Furthermore, observe that 



The proof draws on and combines a number of results from the theory of weak convergence of probability measures 
that can be found in Ethier and Kurtz [15, Chapters 3, 6] 
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= ■.S<U}C {X,(")(a;) : s < Uz}. (1.40) 



Since {X("'(a;)}„gN is convergent in £'r[0,(X)), it follows that (I) for every rational t € [0,U] there exists a 
compact set Ct cM. such that Y^"'\uj) e Ct for all n and (II) for every U > 0, lim^-^o sup„ ■u;^([/) = and 
as a consequence^^ is relatively compact in £']r[0,oo). 



1.4 European and volatility derivatives 

1.4.1 Call and put options 

We first turn to the valuation of a call option price. In the model under consideration a closed form expression 
is available, in terms of the original parameters, for the Fourier transform in log-strike k = log K of the 
call prices Ct{K) with maturity T, 



c*j.{^) = j e'^^CT(e'=)dfc where 9(0 < 0. 



Proposition 7 Define for any ^ e C\{0, i}, a: e R and j e the value D{^, x, j) is defined by the formula 

D{^,x,j) := • [exp{r(if(l + - Ad)]i] (j). (1.41) 

Then if 9(^) < holds that 

where x = log is the log-price at the current time and Zq = j the initial level of the volatility. 

Remarks, (i) The call option price can now be calculated using the method described in Carr-Madan [13] 
by evaluating the integral 

CT{K)='-^^t^ r e-^^'c*As-ia)ds 

^exp(-afc) r ji[e-i«fe^(^_i^^iog5^^^^)]ds, (1.42) 
Jo 

for k = log(_ftr) and any strictly positive a. The integral in (1.42) can be approximated efficiently by a finite 
sum using the FFT algorithm (sec [13]). Since in our model we have an explicit formula for the transform 
c^{s) given by (1-41), the pricing of European call options is immediate. 

(ii) A simple calculation shows that the put option price Pt{K) = \{B!^)~^{K — St)'^\ can be 
expressed in terms of the formula for D{^,x,j) in (1.41) and any strictly negative constant a in the following 
way: 



Both applications follows from the fact that conditions (I) and (II) are necessary and sufficient for the relative 
compactness of (F^"^)„gN (Ethier and Kurtz [15, Theorem 3.6.3]). 
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Pt{K) = ^^P^""'^) f ^\e-'-'''D{s- ia,\ogSo,Zo)]ds, where fc = log(i^). 
Jo 

Proof. To find the European call option price CxiK) = E^^ [(J5|')-1(5t - K)+] in model (1.19) we first 
need to find the Fourier transform in the log-strike k = log(iir) of the function 

CT{k) = exp(afc)E^j [{B^)-\St - exp(fc))+] , 

where a is some strictly positive constant. Fubini's theorem and the form of the characteristic function (1.31) 
imply the following for ^ = v — ia: 



4(6 



[ exp((iw + a)k)E^,j [{B^)-^{St - exp(A;))+] dk 
Js. 

{B^)-^ [ exp{{iv + a)k){ST-exp{k))+dk 
Js. 



= E^j [(Bt exp((l + a + ±v)Xt)] /{a^ + a-v^ + i{2a + l)v) 



px(l+a;+lt>) 



a'^ + a — v'^ + i{2a + l)v 



[exp(T{K(l + a + ±v)- Ad))1] (j). 



This concludes the proof. 



□ 



1.4.2 Implied volatility at extreme strikes 

The implied volatility ax,i{K,T) for a given strike K and maturity T is uniquely defined by the identity 

C^S(5p^i^^T,a,,,(A',r)) =E„ [{B^^)-^(St - K)+] , (1.43) 

where C^^{So, K, T, <t) is the Black-Scholes formula and 5*0 = exp(.T). The results in Lee [23] and refinements 
in Benaim and Friz [6] imply that in model (1.19) the slope of the volatility smile is uniquely determined by 
the quantities , i = 1, . . . ,n, defined in (1.21). In the particular case where the distribution of jumps is 
double exponential, af are in fact the reciprocals of mean-jump sizes in model (1.19). 
In order to state the precise result, define 

q+ := sup {u : E^^i [5^+"] < oo for all i € E°} , 
q- := sup {u : Ex,i [S^"] < oo for all i G E^} . 

If the chain Z is irreducible, the quantities q± can be identified explicitly to be equal to 

q+ = min{af - 1 : i e {1, . . . , iVo} & pAi > 0}, (1.44) 
q_ = min{ar : i e {1, . . . , No} k {1 - pi)Xi > 0}. (1.45) 

As noted above the quantities g+ and g_ depend only on the mean-jump sizes of the compound Poisson 
processes in model (1.19). 

Denote the forward price by Ft := Ej:_i[6'T]. Then the asymptotic behaviour for the implied volatility is 
described as follows: 
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Proposition 8 Suppose that Z is irreducible. For T > and K > and with q± given in (1.44) -(1.45), it 
holds that 

K^oc \0g{K/FT) \y + 



if^o \\ogiK/FT) 

Remark. Note that, if the Markov chain Z is irreducible, the asymptotic slope of the implied volatility smile 
for large and small strikes depends neither on the spot So = e^ nor on the starting volatility regime i. 

1.4.3 Forward starting options and the forward smile 

A forward starting call option is a call option whose strike is fixed at a later date as a proportion of the 
value of the underlying at that moment. More precisely, the pay-off of a Ti-forward starting call option at 
maturity T2 > Ti is given by 

Denote the current value of this forward starting option by Ft^ ,t2 (k) and let 

f¥i,T.(0= /e'«'=FT„T.(e'=)dfc, where $>(0 < 0, 
Jm. 

be its Fourier transform in the forward log-strike k = log k. 
Proposition 9 For ^ with 3(^) < it holds that 

P^MO = ^^-[exp{T,{Q-Ap))exp{{T2-T,){K{l + iO-AD)}l] (j), (1.46) 
where x = log is the log-spot price and Zq = j the initial level of the volatility. 

Remark. An inversion formula, analogous to the one in (1.42), can be used to obtain the value Fti,t2('«) 
from Proposition 9. 

Proof. The price of a Ti-forward starting option is given by the expression 

The process {X, Z) is Markov and therefore, by conditioning on the cr-algebra generated by the process up 
to time Ti, employing the form (1.31) of the characteristic function of Xt and using the spatial homogeneity 
of the log-price Xt = log St in our model, we obtain the following expression for the price of the forward 
starting option 
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(1.47) 



St, 



5' 



Ti 



^^A{B?,-tX\St,-t,-k)+] 



Soe[cxp{T{K{-i) - Ad))Ct,-t,{k;1), 



(1.48) 



where Ct2-Ti {k; 1) is a vector (of call option prices) with j-th component equal to Eoj[(-B^_j.^) ^{St2-Ti ■ 
«)+]. The martingale condition in (1.25) and Proposition 7 conclude the proof. 



□ 



Remarks, (i) A quantity of great interest in the derivatives markets is the forward implied volatility 
^xUSt, k,T) at a future time T implied by the model. It is defined as the unique solution to the equation 



{St, , kSt, ,T2-T,, ai^;, {St,,k,T,))= 



||(5t,-«5tJ + 
T2 



Sti 



(1.49) 



where the left-hand side denotes the Black-Scholes formula with strike kSt, and spot St, ■ The reason for the 
importance of the forward implied volatility (tI^{St, k, T) lies in the problem of hedging of exotic derivatives 
using vanilla options. If at a future time T the spot trades at the level St, then the trader needs to know 
where, according to the model, would the vanilla surface be trading at. This is of particular importance when 
hedging a barrier contract that knocks out at the level St, because conditional on this event the trader is 
left with a portfolio of vanilla options that was created as a semi-static hedge for the exotic derivative. 

(ii) In the model given by (1-18) we can compute the right-hand side of (1.49). In view of the Markov 
property of the process {X, Z), this equation is equivalent to 

CBS(5r^,^c5ri,T2 -ri,afj(^Ti,K,Ti)) = ST,E,,i [Eo,z,^ [{B^.^tXHSt^-t, - k)+] \ St,] 



= St, Y1 ^-'^ [^^1 = ^ 

= ST,r\XT„T^)'CT2-T,{K,l), 



where the coordinates of the vector f^'^{y,T) are defined by 



Zt = j 



XT = y 



(1.50) 



and Ct2-t,{h', 1) is as defined in the line following equation (1.48). 

(iii) The vector Ct2-Ti(k;1) can be computed by formula (1.42) and, in the light of definition (1.50), 
Proposition 3 and formulae (1.29) and (1.30) for q^^{y,j) and q^^{y), it follows that 

rXAf rr>\ xAf \ xAf -\ 

jj {y,T)qT [y) = Qt {y^J}- 

This yields the quantity in (1.50) and hence a formula for the forward implied volatility in our model. 
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1.4.4 Volatility derivatives 



An option on the realized variance is a derivative security that delivers c/j^St) at expiry T. where c6 : R+ — )■ M 
is some measurable payoff function and St is the quadratic variation up to time T of the process logS* = X. 
More formally, for a refining sequence of partitions'^ (-ffn)neN of the interval [0,T], St is given by 



St ■■= lim V log ( ) 

tfenn,i>i \ '-1/ 



It is well-known that the sequence on the right-hand side converges in probability, uniformly on compact 
time intervals (see Jacod & Shiryaev [20], Theorem 4.47) and the limit is given by 

f-T 



St= f a{Z,fdt + E E hz.=i}{Mf- 



(1.51) 



where AJl := .][ — Jl_. The process {St\t>o is called the quadratic variation or realized variance process of 
X, and its law is explicitly characterised as follows: 

Proposition 10 (i) The process {{St, Zt)}t>o is a Markov process with 



St = f a{Z,fAs f I{z.=i}dJl 

Jo .ggo^O 



where , i & E'^ , is a compound Poisson process with intensity Xi and positive jump sizes Ki with probability 
densUy ^ 



2v^ 



^/3+e^^^(-B+)l + (1 -p,)/3-e^^r(_i3-)i I^^^^^{x). 



(a) The discounted Laplace transform of St is given by 

[ ^^P(-^^*) 1 = [e^p{tiKj:{u) - An))l] (i), u>0, 
where Kz;{u) = Q + As{u) with As{u) an Nq x Nq diagonal matrix with i-th element given by 

V'f (u) := -uaf + Xi (E [exp(-ui^i)] - 1) , 

with 

E [exp(-«i^,)] = {^^^) + ~ ^'^^'^ {7^^') ^"^'"0 ^' 

where <P{x) := exp(a;^/2)A/'(a;), with the cumulative normal distribution function J\f . 



(1.52) 



(1.53) 



The sets i7„ = {to , i", . . . , tJJ }, n e N, consist of increasing sequences of times such that to =0,tn= T, n„ C i7„+i 
for all n e and lim„^oo max{|t" — | : i = 1, . . . , n} = 0. 
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Remarks, (i) As a given matrix M in practice typically^^ admits a spectral decomposition M = UDU~^ 
where Z) is a diagonal matrix, <?(M) can be evaluated by ^(M) = U^{D)U~^ where <?(-D) is the diagonal 
matrix with i-th element <P{Dii). 

(ii) It is important to note that the realized variance process S does not possess exponential moments of 
any order. This follows directly from the fact that the distribution of jumps gi given in (1.60) decays at the 
rate e~'^^, for some positive constant c, and implies that the left-hand side in formula (1.52) will be infinite 
for complex numbers u with negative imaginary part. 

(iii) The expression (1.52) for the discounted Laplace transform can be employed to obtain explicit results 
for the values of volatility derivatives. The buyer of a swap on the realized variance pays premiums at a 
certain rate (the swap rate) to receive at maturity a pay-off 4'{St) that is a function </> of the realized variance 
Et, with as most common examples the volatility and the variance swap. In the case of a variance swap this 
function is linear {4>{x) — x/T) whereas for a volatility swap it is a square root ((/>(a;) = \/xJT). The swap 
rates are determined such that at initiation the value of the swap is nil. 

Corollary 1 Suppose that Zq = j. Then the variance and volatility swap rates ^var{T,j) and ^voi{T,j) are 
given as follows: 

W(T, j) = r { [exp(T(Q - .1^)) - exp{T{Kj:{u) - An))] l}{j)^, (1-54) 

w(T,j) = [exp{t{Q - AD)}AveM{T-t){Q - ylc)}l](j)di (1.55) 

= ^ [{exp(T(Q - An)) - eMT{Kj:(h) - An))} 1] (i) + o{h), /i 4 0, (1.56) 
where Ay is a Nq x Nq diagonal matrix with i-th element given by 

V{i) = + 2A, (p,(^+)'(B+)-2 + (1 - pO(/3r)'(Sr)-') 1- (1-57) 
Remarks, (i) The Laplace transform ayariQ^j) of cr^ar(') j) • T T<;yar{T,j) is explicitly given by 

^var{q,j) = [{ql -An- Q)-'Av{qI -An- Q)-'l\ {j). (1.58) 

(ii) It is clear from the definition of Ks{u) that the integral in (1.54) converges at the rate proportional 
to 1/VU where U is an arbitrary upper bound used in the numerical integration in (1.54). 

Proof of Proposition 10. It is clear from the representation (1.51), that the increment Ef — Sg, for any 

t > s >0, satisfies the equation 

St-Ss= f <j\Zu)du +Y. I' hzu^i}'^Ju. (1-59) 



This is the case since the set of all square matrices that do not possess a diagonal decomposition is of co-dimension 
one in the space of all square matrices and therefore has Lebesgue measure zero. 
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where J^, i G E^, is a compound Poisson process with intensity Aj and positive jump sizes Ki distributed as 
(J/*)^ where Ui follows a DPH{pi, , Bf , , B^) distribution. In particular, Ki is distributed according 
to the density 

1 



9i{x) 



2v^ 



(1.60) 



where the probability density function fi is given by 

n{x)=Pi{^t)'^''''H-Bt)'^I(o,oo){x) + {l-Pim)'e-''^-Br^^^^ 

As the J* have independent increments, and Z is a Markov chain, it directly follows from (1.59) that {St, Zf) 
is a Markov process, and moreover, a regime-switching subordinator. The form of the discounted Laplace 
transform can be derived as in Theorem 1. ^ 



Proof of Corollary 1. Employing the following representation for the square root 

1 



20F Jo 



1 — exp(— ua;)] 



3/2' 



for any a; > 0, 



as well as the form of the discounted characteristic function given in (1.52) and Fubini's theorem yields that 
the volatility swap rate can be calculated via a single one-dimensional integral 



B^ 



1 



2VttT 



i [exp(T(g - Ad)) - exp{T{Ks{u) - Ad))] 1 



du 



The derivation of the variance swap rate formula (1.55) rests on a conditioning argument. Indeed, by condi- 
tioning on the sigma algebra = <j{{Zt}t<T) generated by Z up to time T, it follows that 



E, 



{Jo .^„oJo 



1 
Jo 



where w{j) := a'^{j) From the definition of it is easily checked that E[Jf ] is equal to Xj times the 
second moment of the density /j. One verifies by a straightforward calculation that w(j) is equal to V{j) 
given in (1.57). Furthermore, the Markov property of Z applied at time t yields that 



E.. 



-Dt Jo 



e'i exp(t(Q - AD))e',ej exp((r - t){Q - Ad))1. 



Equation (1.55) follows then by an application of Fubini's theorem. Furthermore, Equation (1.56) follows by 
noting that the expectation ^x,i[^T/B^] can also be obtained by calculating the negative of the derivative 
of the Laplace transform at zero: 



--{[eMTiK^iu)-AD))l] (i)} 



u=0 



□ 
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1.5 First passage times for regime-switching processes 



1.5.1 Three key matrices 

The characteristics of the process {X, Z) can be summarised in terms of three matrices Qq, S and V that 
will shortly be specified. Given those three matrices we will show how to reconstruct {X,Z) in Section 1.7. 

The matrices Qo, ^ and V will be specified in nine- block matrices using block- notation; the middle block 
of Qoi ^ and V describes the rates of regime-switches, and the volatility and drift of the process in the 
different regimes, while the upper left and lower right block of Qq specify the distribution of up-ward and 
down-ward jumps in the different regimes, in terms of the (phase-type) generators. More precisely, we define 
the three key matrices Qq, S and V, in block notation, by 

( B+ b+ 0\ 

A+ Q-Ax A- , (1.61) 
b- B- 



S := 





Qo ■■= 


(o 


O 





As 


\0 


O 



o 

o\ 
o 

0/ 



I 



V := 








o\ 


o 


Am 


o 









Am 


denote A^o 


X Nq 



(1.62) 



elements 



Ax{i,i) := \i, As{i,i) := u{i) and AMii,i) ■= Ki)- 



(1.63) 



Further, O and / are zero and identity matrices of appropriate sizes such that Qo, S and V are square 
matrices of the same dimension. In block-notation A^ , and are given by 



A^ := 



\ 



\ 



\ 



\ 



6± := 



(-Bfl 



\ 



\ 



(1.64) 



where := XiPi and := Ai(l - Pi). 

Remark. The matrix Qo is in fact the generator matrix of a Markov chain, as it has non-negative off- 
diagonal elements and zero row sums. We denote the state-space of this Markov chain by E. In the sequel 
we will frequently use the following partition of the set E: 



E^ 



{i€E:Eiii^Q}, E+ = {i€E:Sii=0,Vii>Q}, E- = {i € E : Su = Q,Vu < Q). (1.65) 



Note further that E^ can and will be identified with the state-space E^ of the chain Z. See Section 1.7 for 
further properties of the Markov chain defined by the generator Qq. 
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1.5.2 Matrix Wiener-Hopf factorisation 



For a given vector of discount rates h : E ^ C, the matrix Wiener-Hopf factorisation associates to the 
matrix 



where is a diagonal matrix with i-th diagonal element Ah{i,i) = h(i), a quadruple of matrices which, as 
we will show below, characterises the distributions of the running maximum and minimum of X. 

Let us briefly describe the sets of matrices of which this quadruple are elements. Denote by P(n) the set 
of n X n square matrices whose eigenvalues all have non-positive real part. Note that by Lemma 2, P(n) 
includes the set G(n) of n x n sub-generator matrices (i.e. matrices with non- negative ofF-diagonal elements 
and non-positive rows). Recall that C"^™ denotes the set of n x m matrices with complex entries. Denote by 
IHI(n, m) the set of n x m sub-probability matrices (i.e. matrices with non-negative elements and row sums 
smaller or equal to one). 

Denote by A^, N+, N+, N' and N~ the number of elements of the sets E, E^ U E+ , E+ , E° U E' , and 
E~ , respectively. Also, let H denote the set 



Definition. Let he V. and let W+,G+, W' and be elements of the sets C^x^"", D(7V+), C^xJV 
©(TV^), respectively. A quadruple (W^+, G+, ,G^) is called a matrix Wiener-Hopf factorisation of Qh if 
the following matrix equations are satisfied: 



where 0+ and O are zero matrices of size N x and N x N . 

Theorem 4 (i) For any h £H, there exists a unique matrix Wiener-Hopf factorisation of Qh, denoted by 



(a) If h ^ ^R{h), then G G(iV+) and G G(iV~) are sub-generator matrices, and rj'^ G H(iV, iV+) 
and rj^ e W{N,N~) are in block-notation given by 



Qh ■= Qo - Ah 




min mh(i)) > 0, min mh(i)) > 




(1.67) 



(1.66) 




and 




(1.68) 



for som,e mMrices e M.{N_ , A''+) and r]^ G H(7V"'", N ), and identity matrices 7+ and 7_ of sizes x 7V+ 
and N- X N-. 



Remarks, (i) For h given by h{i) = qlE°{i) with g > we will also write {Vjq ,Qq ,r]q ,Qg ). The proof of 
Theorem 4 will be given in Section 1.8. 
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(ii) We allow the vector h to take complex values to be able to deal with a Laplace inversion using a 
Browmich integral, which involves the integration of the resulting first-passage quantities over a curve in the 
complex plane. Note that, for any real- valued h G H, i.e. h = the matrix Qh is a transient generator 

matrix, since the off-diagonal elements of Qh are non-negative and Qftl ^ 0. 



1.5.3 First-passage into a half-line 

The marginal distributions of the maximum and minimum as well as the distributions of the first-passage 
times into a half-line can be described explicitly in terms of the matrix Wicncr-Hopf factorisation. Denote 
by Xf = supo<s<j Xs the running maximum of X at time t and by X.t — info<s<t Xs the corresponding 
running minimum, and let and T~ be the first passage times of X into a half-line, 

T+ = M{t > : Xt e (o, oo)}, T" = M{t > : Xt £ (-oo, a)}. 

The distributions of those random variables are related via 

r+it) = P.,,(T+ < t) - P,,,(Xt > a), 

r-{t) = p.,^(^- <t) = F^4-Xt > a). 

If Bq denotes a random time that is exponentially distributed with parameter q > and that is independent 
of {X, Z) , then the running maximum and minimum at time follows a phase- type distribution, with 
parameters explicitly given in terms of the matrix Wiener- Hopf factorisation (77+, Qq,ri~, Qg). 
For any h gH, define the N x N+ and N x N~ matrices ^^^^(x) and ^^a^^) ^J" 

<a(a^'^j) = Kexp((a-a:)Q+)] (^,JV(-oo,a](a;)+5^J/(a,oo)(a;), i G E, j G E+ U E"" , (f.69) 
^h,ai^'i'j)= feexp((a;-a)Q;;;)] (i,j)7[a,oo)(a;)+5ij/(-cx.,a)(a;), i G E,j G E' [J E° , (1.70) 

where dij is a Kronecker delta (i.e. Sij = I^iy{j)). For h given by h{i) = qlEo{i) with g > 0, we will also 
denote ^h,ai^) by ^g,aix)- 

Proposition 11 For q> it holds that under F^^i 

X^^-x^ PH{ij+{i),Q+), xGR,iGE°, 
-2L^+xr^ PH{n-(i),Q-), xGR,iGE°, 

where r]^{i) and r]~{i) are the i-th rows of r]g and r}~ . In particular the Laplace transforms r^{q) ■= 
e~'^*r'^{t)dt, for q> 0, are given by the formulae 



^^(9) = ^ E ^^A^-'^' hz^^=j}] = ^ E '^U^'^'j)' x,aGR,iG E\ 
^"(9) = - E Ex,i[e-''^»"/{z^_=,}] = - 5Z ^^a(a;.«,i). x,aGR,iG E°. 



q — ' ^ ' q 
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Remarks, (i) The proof of Proposition 11 will be given in Section 1.8. 

(ii) Denote by i) and ^~^^{x,i) the A''+-dimensional and A''" -dimensional row vectors with j-th 

elements <?+„(a;,i,j) and <P~ ^{x,i,j), respectively. Then Proposition 11 implies that, under the probability 
measure Px,i, x gM., i & the processes M+ = {M^}t>o and M~ = {M^}t>o, defined by 

are row- vectors of bounded martingales. Indeed, the Markov property of {X, Z) implies that 



= I 



iX,Z) 



= e " 'i!^x +,z + 



= e+(i) if 



where denotes the j-th standard basis vector in and we used that K^^i e~'^^-" I. 
X > V and i G i?*^, which directly follows from the definitions (1.69) and . 

(iii) Suppose that Qo = O. This corresponds to a model in which there are no jumps and no switches 
between the regimes (i.e. with probability one the process stays in the starting regime and evolves as a 
Brownian motion with drift). In this case we can identify the matrix Wiener-Hopf factorisation in closed 
form. Note that we have = N~ = N and hence the matrices and are of dimension N x N.lf we 
take to be equal to the identity matrix and h{i) = Roii) + Q, the matrix equations (1.66)-(1.67) reduce 
to 

IjJ^G+f - VG+ = {Ad + ql) = \s\G-f + VG", 

where / is an iV x TV identity matrix (recall that the discount rate matrix Ajj is diagonal and satisfies 
Aoih i) = Roii) for all i e E^). These equations are satisfied by the diagonal matrices 



G+ = diag(-a;,+ , i = l,...,N), G" = diag(-a;. , i = 1, . . . , iV), 



(1.71) 



where 



± Ml 



Ml 



2(<Z + ri) 



and Ti := Roii)- 



(1.72) 



In particular, we obtain the well known fact that the maximum of a Brownian motion with drift at an 
independent exponential time is exponentially distributed: 



(Xe,+,^ > a) = e-"^(»-") for x < a. 



32 Aleksandar Mijatovic and Martijn Pistorius 

1.5.4 Joint distribution of the meLximum and minimum 

In the previous section we have shown how the marginal distribution of the maximum as well as of the mini- 
mum can be explicitly expressed in terms of the matrix Wiener- Hopf factorisation. Also the joint distribution 
of the running maximum and minimum, 

can be explicitly identified in terms of the matrix Wiener-Hopf factorisation, by considering appropriate 
linear combinations of the functions and defined in (1.69). 

To formulate the result, introduce the matrices Z+ e and Z~ e C^^^^ by 



i€E°yjE+, jGE^UE-, 



and define, for any h gH and any a; € M, the N x and N x N matrices ^{x) and ^ ^{x) 

^M,.(^) = (r/,;c«^"(^-'' - ry+eQ^(---)Z-) (/ - Z+Z-y' 
where A'^ and A~ are N x and N x A''" matrices with elements 



(1.73) 
(1.74) 



A+{i,j) = 5ij, A {i,k) = 6ik, 



z e i e u fc e £° u E- 



For h given by h{i) = qlEo{i) with g' > 0, we will also write ^{x) = ^^(^,ix). 
Proposition 12 Let i G E^ and 

^liit) = P,,i{Ti,^ < t, Xt,,„ > v), i^-^S) = ^xATt,v < t, Xt,,„ < £), 

where 

T^,„ := inf{i >0:Xt^[i,v]} = Tf A T+ . 
For any q > the Laplace transforms in t of^^l)'^^ and are given by 

The Laplace transform of tpx,i{t) = Px,i{X_t > (-^Xt < v) is hence of the form 

V^x,i(g) =V^+i(g) + C,i(9)• 



(1.75) 
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Remarks, (i) If Qo = O, then r}^ are identity matrices and the following identities hold: 

= (gQjC^'-a;) _ QQ-{x-e)^Q+{v-l)^^ _ ^Q-{v-t)^Q+{v-t'.)yl^ 



where are diagonal matrices given in (1.71). In particular, we find the well known two-sided exit identity 
for Brownian motion with drift: 

where ujf^ is given in (1.72) with ri ~ 0. 

(ii) In the case that no volatility is present {U = 0) the identities simplify and we find the expressions 
in [5]. 

(iii) Under P^^^^, x e M, i € E° , the process M+ = {Mj^}t>Q defined by 

M+ = e-«(*^^^-)!Z'+ (X,^r,,„,^tAT,,J , 

is a row- vector of bounded martingales, where we denoted by ^qi^{x,i) the A''+-dimensional row vector 
with j-th element 'l'q£^{x, 

Proof. Let g > and define 

In view of the fact that M+ is a bounded martingale, it holds that 

g+{x,i)= Yl ^-,i[^~'"K,i,vi^r,Zr,j)I{r<oo} 

= ^^'i [e~''"I{x^>v,z^=j}] 

= q / e-«'P,,.(r <t,Xr> v)dt = q^^tM- 
Jo 

Here we used that, in view of the definitions (1.73) and (1.68) of ^ and the function 5+ satisfies: 

g+{x,i) = l i{x>v,iGE° 
g+{x,i)=0 iix<e,ieE°. 

The expression for can be derived by a similar reasoning. ^ 
1.5.5 Valuing a double-barrier rebate option 

A double-barrier rebate option pays a constant rebate L at the moment r one of the barrier levels is crossed, 
if this happens before maturity T. By standard arbitrage pricing arguments the Laplace transform Vreb in 
maturity of the price of Vreb such an option is given by 
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'cxp(-rjr) 



We will find below the following more general quantity, which will also be employed in the sequel: 



H^'\q, u, j) := E,,i \{B^) ' exp {±uX^ - qr) I^z^^jy 



for j = l,...,N. 



(1.76) 



(1.77) 



In the following result an explicit expression is given for the quantity J'* in terms of the matrix Wiener-Hopf 
factorisation, which is an extension of the identity (1.75). 

We denote by and , for i e E''' = {1, . . . , N}, the parts of the state-space E corresponding to the 
blocks and in the matrix Qo in (1.61). Recall that the definition of af was given in (1.21). 

Theorem 5 For any h gH, i,j G E^, and u gC that satisfies ?s{u) G (— ai",aT) it holds that 

Hj'^h, u) = [^Lkl^) (*) + {^u,.klu) (0, (1-78) 
where the column vectors fc^^ = {kj'u{i),i G E^ U E^) and = {k~^{i), i G E^ U E~) are given by 



{{-ul+ - B+)-\-B+)l) ii), 



tfi = jGEO, 

if 7. e E+. 
otherwise 



{{uI.-BJ)-\-B-)l) (i), 
0, 



ifi = j e E°, 

ifie EJ, 
otherwise 



where and Ij are \Ej'\ x \Ej'\ and \Ej \ x \Ej \ identity matrices. 
The proof is given in Section 1.8. 



1.6 Double- no-touch and other bsirrier options 

In this section we show how the prices of doublc-no-touch and double knock-out call options can be expressed 
in terms characteristic function H of the process at the corresponding first passage time, which was identified 
in Theorem 5. 

A double-no-touch is a derivative security that pays one unit of the underlying asset at expiry T if the 
underlying asset price does not leave the interval [exp(^),exp(i;)] during the time period [0, T], where £ <v. 
Similarly a double knock-out call option struck at K delivers the payoff {St — K)~^ ■= maxjS'T — K, 0} if 
throughout the life of the option the asset price stays within the interval [cxp{£),cxp{v)]. The arbitrage free 
prices for a double-no-touch and a double knock-out call option are respectively given by 

'hr>T}' 



Cx,i{k, T) = 



B^ 

hr>T} 



B, 



(1.79) 
(1.80) 
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where k = logK is the log-strike and r is the first time the process S leaves the interval [exp(^),exp(i;)] or 
equivalently 



T := mi{t > : Xt [i,v]}. 



(1.81) 



We will find it more convenient to consider the double-touch-in and the knock-in call whose values Vdu and 
Vkic, as function of maturity T and log strike k = logi^^, are given by 



Vdti{T) := E^,i 
Vkic{T,k) := E^^i 



'{t<T} 



(Si?) \sT-e'')+\-C,,i{k,T). 



(1.82) 
(1.83) 



Since the zero-coupon bond price as well as the call option price have already been identified, the problem of 
calculating the double- no-touch and knock-out call prices thus reduces to identifying Vdu and Vkic- In general, 
no closed form expressions are known for these function in terms of elementary functions. Below we show 
that the Laplace transform in T of Vdti{T) as well as the joint Fourier-Laplace transform in {k,T) of Vkic 
can be identified explicitly in terms of the parameters that define the log-price process X . Both transforms 
will be identified in terms of Laplace transform Fx.i{u, q) in T of the function T t-^ Fx^i{u, T) given by 



F^^i{u,T) :=E^,i 



qd exp(iuXT) 



(1.84) 



The Laplace transform Fx^i{u, q) in turn will be given in terms of the A^o-vector H^'^{q, u) whose coordinates 
are given by 



m'\q,u) := E, 



e^p\^±uXr-J {RoiZs) + q)dsj I{z^=j} 
Theorem 6 For any q> and ^ with S>(^) <Q it holds that 



for j = l,...,Aro. 



vdtiiq) = -P'x,i(o,g) 
1 



F,4^-±,q). 



(1.85) 

(1.86) 



Here, Fx,i{u,q) is given by 

F,4u, q) = {H^^\q, u))' {ql + Ad - K{u))-H, (1.87) 
for all q &C, such that 5R(g) > q* = max{5i(V'i(w)) - i?r>(«) : i = l,..., Nq}. 

Remarks, (i) Note that if u = 0, then q* < 0, as the interest rates Roii) is assumed to be non-negative, 
(ii) The Laplace transform in (1.87) can be inverted by evaluating the Bromwich integral 



-I pc+±oo 

F,,i{u,t) = — / e"i{H-'\q,u)y{qI + An-K{u))-'ldq, 

^'^ Jc— too 



(1.88) 



36 Aleksandar Mijatovic and Martijn Pistorius 

for any c> q*. An efficient algorithm to approximate this integral can e.g. be found in Abate and Whitt [1] . 

Proof. First note that by a similar calculation to the one in the proof of Proposition 7 we find that the 
Fourier-transform in log-strike k of e"''vkic{T, k) can be expressed in terms of F by 



/ 



Vkic{T, k)dk 



1 



(a + iv){l + a + iv) 



for any a > 0. 

Assume next that g > is real with q > q*. The Laplace transform Fx^i{u, q) has the following well-known 
equivalent probabilistic representation 



(1.89) 



where Cg is an exponential random variable with parameter q which is independent of the process {X, Z). 

By conditioning on J^t and applying strong Markov property at r of the process {X, Z) together with 
Lemma 4 we find that 



= IE.. 



exp(iMX^) 
B? 

exp(i'uXT-) 



E 



-f{r<e,} exp 



{Z,)&S + iu(Xe, - Xr) 



^Hj'\q,u)h{j,q,u) 



h{Zr,q,u)e 



(1.90) 



where the value h{j, q, u) for each state j G of the Markov chain Z is given by 



Kj, q, u) := Eo„ 



exp 



RDiZs)ds + iuXe 



[{ql + Ad - K{u))-'l] (j) 



(1.91) 



for some exponential random variable e'g with parameter q that is independent of the Markov process {X, Z), 
where the second equality follows from Lemma 3 and Theorem 1. 

A key observation that follows from the representation (1.90) is that the function 

q^E[F^^i{u,eq)]/q 

has a holomorphic extension to the complex half-plane {(j G C : 3?((j) > q*} which therefore^® coincides on 
this domain with the Laplace transform. Thus, Equation (1.90) holds for q in this domain, and the proof is 
complete. ^ 



If two holomorphic functions defined on a connected open set in C coincide on a subset with at least one 
accumulation point in i?, then they coincide on the entire f2. For a proof of this well-known statement see [29], 
page 208, Theorem 10.18. 
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In this section and the next we wih provide a proof of the matrix Wiener-Hopf factorisation and its corollaries 
derived in previous sections. We will proceed in two steps: 

(i) Reduction of the first-passage problems of X over constant levels, which will involve overshoots and 
undershoots due to the jiimps of X . to the first-hitting problem of a constant level by a regime-switching 
Brownian motion, employing a classical argument for the embedding of phase-type jumps (see Asmussen 
[2]), which we review in this section. 

(ii) Solution of the first-passage problem of regime-switching Brownian motion via a characterisation of the 
dynamics of the ladder processes, which is carried out in Section 1.8. 

Let y be a continuous time Markov chain with finite state space EU d and generator restricted to E 
given by Qo, where d is an absorbing cemetery state and E and Qo will be specified shortly, and denote by 
^ = {^t,t > 0} the Markov modulated Brownian motion given by 



where a; e IR is the starting point of ^ and s and m are functions from E U d to M. also to be specified 
shortly. The couple (^, Y) defined as such is a two-dimensional strong Markov process. In the sequel we will 
denote by F^.i and 'Ex,i the conditional probability P^.i = lP[-|^o = x,Yq = i] and conditional expectation 
'Ex,i = E[-|^o — x,Yo = i] respectively. 

Let the state-space E be as in Section 1.5. In other words E = E~ U i?^ U E^ where E^ is the state-space 
of the chain Z and £J+ and E~ are given in (1.65). The generator Qo is given in (1.61) and m{i) := Av{i, i), 
s{i) := Aj:{i,i) for all i G E, where the matrices Ay, As are defined in (1.62). Thus, while the chain Y 
is in state j e E^, ^ evolves as a Brownian motion with drift m{j) and volatility s{j), and while Y takes 
values in and E~ , the path of ^ is linear with slope -|-1 or —1. Informally, a path of X can be obtained 
from a path of ^ by replacing these stretches of unit slope by jumps of the same length, as is illustrated in 
Figure 1.1. These linear increasing and decreasing stretches of path of ^ thus correspond to the positive and 
negative jumps of X, respectively. 

More formally, by appropriately time-changing (^, Y) a stochastic process can be constructed that has 
the same law as {X, Z). Denote by 



the time before t spent by the chain Y in E^ and its right-continuous inverse respectively. It is clear from 
the definition of the generator Qo in (1-61) that when the chain Y jumps from any of the states in E^ U E~ 
to a state j G E^, it must have been in the state j just before it had left E^. Therefore the form of the 
matrix (1.61) implies that the process Y o Tq^, which simply ignores all excursions of Y into E~^ U E~ , is 
a Markov chain with the same generator as Z. Furthermore it is straightforward to verify that ^ o Tq^ is a 




(1.92) 
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Fig. 1.1. Shown is a sample patii of X until the first time t that X exits the interval [l,v\. The process ^ has no 
positive jumps and always hits a level at first-passage. 



regime-switching jump-diffusion and the law of the process (C ° Tq ^ ^Y °Tq ^) under P^^^ is equal to that of 
(X, Z) under P^,i, for all x e M and i G . 
If we define the first-passage time 

T« = inf{t > : ^ v] and Yt G E^} (1.93) 

it follows that random variables ro(r^) and the stopping time r defined in (1.81) have the same distribution. 
In particular, with the extension of h to E that puts h{i) —{) ioi i ^ E^ and that we will also denote by h, 
it holds that 

( ^Ti^ / /i(l^s)ds, l^e I under P^.i has the same distribution as ( Xt, / h{Zs)As, Zr ) under 



for X e M and i G E^ . The function iJj'* defined in (1.77) can thus be expressed in terms of the embedding 
(^, y) as follows: 



exp ( - 1 h{Ys)(ls I /{y^j=j} 



i,je E^. 



(1.94) 



1.8 Ladder processes 

This section is devoted to the proof of existence and uniqueness of the /vmatrix Wiener-Hopf factorisation, 
for any h H. For real h £ H the vector h has the probabilistic interpretation of a vector of state-dependent 



1 Exotic derivatives under stochastic volatility models with jumps 39 

rates of discounting, and in this case the matrix Wiener-Hopf factorisation has a probabiUstic interpretation 
in terms of ladder processes. The existence and uniqueness results are extended to the case of general complex 
h hy analytical continuation arguments. It is of interest to consider the case of complex entries, as the 
Laplace transforms of barrier option prices are expressed in terms of the matrix Wiener-Hopf factorisation 
as we saw above and some widely used Laplace transform inversion algorithms are based on the Bromwich 
(complex) integral representation. 

A classical probabilistic approach to characterisation of the joint distribution of r+ and the position of 
Y at r+ is to consider the up-crossing ladder process of (^, Y), defined as follows: 

Definition. The up-crossing ladder process F+ = {y^"'"}t>o of {^,Y) is given by 

fr(T+), ifT+<00, 

Y+ -.= 1 ^ " (L95) 

I d, otherwise, 

where 9 is a graveyard state and 

r+ = inf{s >0:^s>a} 

with inf = DO. 

Remark. In the case that the original chain is killed at (state-dependent, real- valued) rate h, the up-crossing 
ladder process can be defined as follows. Recall from Markov chain theory that the chain Y^ with state-space 
EU {d} and generator matrix 

('o'* o'* ) where Qh '■= Qo — Ah and is the diagonal matrix with {Ah)n = hi 

has the same distribution as the chain Y killed (i.e. sent to the graveyard state d) independently at rate h{i) 
when Yt = i. In particular, for i,jGE and t>0 it holds that 



„-/'/j(K,)dsr 

e Jo y > ^{Yt=j} 



P(r,^ = j\Yo^ = i)= eMQhtKhj). (1.96) 



If we denote by ^'^ = {^t}t>o the process defined by (1.92) driven by the "killed" chain Y^ instead of Y, 
then the up-crossing ladder process Y^~^ is given by 



Y\+ if tI^+ < DO, 
d otherwise, 



where r^'"*" is the stopping time defined by 

= inf{s >0:C^>a} 

with inf = 00. 



Proposition 13 The process y+ is a Markov chain with state-space E° U E^ U {d}. 
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Remarks, (i) We denote by Qq the generator of F+ restricted to E^[JE~^ and by tjq" the initial distribution, 
given by 



71+ ii,j)= Po.i [^0+ = J, To+ < oo] for ieE,jeE+UE° 



(1.97) 



(ii) More generally, for any real valued ft, G "H it holds that Y'^+ is a Markov chain, with generator, 
restricted to U E+ denoted by and initial (sub-)probability distribution by rj^ = {r]'^{i,j),i & E,j £ 
E+UE°) with 



^ ^ {Y + =j, r + <oo} 



for ieE,je E+UE^. 



(1.98) 



We have the following identities analogous to (1.96): 



E 



{il+=J.n"*'<oo} 



Yo = i 



E 



e I, 



Y, 



{Y^+=j^r+«x>} 

P = i,r/^+< 00 1 F\+=t) 

F{Y/^+=j\Y,^+=i) 



+ =1 



(1.99) 



where we wrote P[-|^] := ¥[-\A f] {Co = 0}] and E[-|^] K[-\A f] {^q = 0}] to simplify the notation. Note 
that the first equality holds since Tq'" = if the chain Y starts at i G E° U E+. In particular. From (1.99) we 
find that fov i,j G E'^ [J E+ 



)+(z,z) =lim-<^Eo,. 



44,0 t 



a. 



(1.100) 
(1.101) 



Proof. Let / be any bounded real- valued Borel function and let b < a. Denote by = (^{Yu}u<a and 
Gt = cr{Ys}s<t the sigma algebras generated by Y+ up to time a and by Yg up to time t. Observing that 
C Grt and using the spatial homogeneity and continuity of ^ and the strong Markov property of Y) 
we find that 



E,,y + [/(>;+)] I 



where we wrote Ej = Eqj. 



□ 
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For any h G H, we define the square-matrix (resolvent) functions u i-)- -R^(w) of dimension on the 
complex half-plane C>o by 



^ ^{r+=j,r+«x)}'^y 
The matrix can then be defined for any h G H, as follows: 
Lemma 5 Let h GT-L and a > and define the matrix G D(Af+) by 



for i,j gE+UE" and u G C>o. (1.102) 



where qi = —Qo{i,i), hi = h{i), Q = Qo + diagjgi : i G E} and 

' fii qi + h 



ifieE+, 



(1.103) 



F 
qi + hi 



ifi G E^ 
ifi G E^ 



where 



Then it holds that 



exp(Q+a)(i,j) i,jGE°UE+. 



{Y+=j,T+<oo} 

Remarks, (i) Note that Proposition 11 follows as a direct consequence of Equation (1.106) since 

P^,i (Xe, >a)= P^,i > a) , 

where ( is the life time of the killed up-crossing ladder process. 

(ii) Equation (1.106) yields in particular the joint Laplace transform of the vector 



(1.104) 



(1.105) 



(1.106) 



Z+ 



I{Y,=i}ds, iGE], 



(1.107) 



whose components record the length of time spent by Y in each of the states in E, until the moment t+ of 

first-passage. 

(iii) The down-crossing ladder process Y~ , defined as the up-crossing ladder process of (— ^,y), is a 
Markov chain with generator restricted to Eo[jE~ denoted by Qq . Analogously to equation (1.103), for any 
h gH a, matrix Qf^ can be defined, which satisfies 



-So" HyMbj 



V^a =3,Ta <Oo} 



exp(g-o)(i,j) i,jGE°UE-. 



(1.108) 
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1.8.1 Proof of Lemma 5 



We will show that the limits in (I.IOO)-(I.IOI) in fact exist for any h GH, and identify these. Let 

p = mi{t >0:Yt^Yo} 

be the first time that the chain Y jumps and — inf{t > : XI > a} the first-passage time of XI := 
Hit + (TjWt (recall that if i e E'^, then /Xj = ±1 and C7j = 0) over the level a and let e, be an exponential 
random time with mean 1/qi (with = —Qo{i,i)) that is independent of X*. In view of the definition of ^ 
(i.e. the first jump time p of F is independent of X'), it follows that 



^{Y+=3,T+<oo,T+<p} 



-Tlh(i)j 



{r;<e,} 



^-Tl(h,+qi) 



= I{i=j,ieE+} exp {-t{hi + Qi)) 

+ -f{i=j,ieso} exp ^-io-r^ + y^/i? + 2(gj + /i^V?^ ^ 

= [1 - ujft + o{t)] as U 0, 



where ujf was defined in Lemma 5 and we used the fact that the Laplace transform of is given by 

= exp ( -aar2 ( + ^fj,j + 2qaA J a>0,qe C>o. 



For every h £ H we can define the following matrices 



® -'{F+=i,r+<oo} 



for i,j e E+UE^. 



(1.109) 



Note that definiton in (1.102) implies that the identity Rl{u){i,j) = e-'^ypl^+{i,j)Ay holds for all 
u e C>o. For every h gH, x & {—oo,t\ and m £ E we have the foUwing identity 



® ^{Y + =i,T+<oo} 



-/o*-»h(n)d.j 

*^ -'{nt.=j.nt.<oo} 



= K j) forieE+U^o, (1.110) 



where jj^ is defined in (1.98). The first equality in (1.110) is a consequence of the spacial homogeneity of the 
process ^. 

Let ^„ = supg<„ denote the running supremum of the stochastic process ^. The strong Markov property 
applied at the first jump time p of the chain Y and (1.110) imply that 



Eo, 



{yt*=J,r+<oo,T+>p} 



meE\{i} 
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The last equality is a consequence of the fact that Yp is independent of the random vector (p, ^p, and takes 
values in the set E\{i} with Fo,i{Yp = m) = Qo{i,m)/qi. Since p is the first jump time of the chain Y, the 
vector (yO, Cpj'fp) has the same distribution as the vector (ej , X* . , . ) , where as above ej is an exponential 
random variable with mean l/qi, independent of the Brownian motion with drift = iiit + aiWt- The 
symbol denotes the maximum of X'^ at the independent exponential time e^. Note that if i € £'+, then 
X^ = XI =t for all t e R+ and the expectation in (1.111) is very easy to compute. 

Assume now that i G E^. Then the Wiener-Hopf factorisation implies that the random variables X^^ 
and 'X^eu ~ exponentially distributed with parameters 

f^i"' (-Mi + V^Mi +2wa2^ and (^fn + ^ ti^^ + 2uaf^ 

respectively and independent for any exponential random variable e„ with parameter u > 0, which is 
independent of X'. Therefore the joint density /-^i y'-x^ satisfies the following identity for any x,y G (0, oo) 

/■OO poo POO 

/ ue-'''UiYi_^,{x,y)dt= / «e-"*/x,Wdi / ue-''' Ui _^,{y)dt for all u>0, (1.112) 

Jo *' t t Jq Jq t t 

where f^i_-^i and /y< are the densities of the corresponding random variables. It is clear that there exists a 
unique extension to the complex half-plane C>o of both sides of the formula in (1.112) and that the following 
formulae hold 



Jo 



Jo 



e-"*/x, (^)dt = — ^ — ^e""''' V-'^'^v ^--^ ar > 0, u e C>o,(1.113) 

ut, . .„,^^._^_V g-y-, (^,+V^?+2-f) y>o,«eC>o. (1.114) 



10 

By substituting the identity (1.112) into the expectation (1.111), applying the formulae in (1.113) and (1.114) 
and taking the limit as t tends to zero we see that the limits in (1.100) and (1.101) are as stated for any 
h gH. Furthermore the formula in (1.103) holds. 

The strong Markov property of {^,Y) next implies that for any h G H and any t > the matrices 
= {P^+{i,j)t, i,j G E^[JE+) with P^+{i,j) defined by (1.109) satisfy the system of ordinary differential 
equations 

^ph+^pH+Q+^ t>0, P^+=I, 

where I denotes an x iV+ identity matrix, the unique solution of which is given by 

= eMQtt). 

This proves that the matrix identified above satisfies (1.106). Since Pt+{i,j) for all i,j G U E+ 
as t — > 00, it follows that all eigenvalues of must have non-positive real part and therefore G D(iV+). 
The proof of the existence of a matrix Qj^ e IS>{N^) satisfying (1.108) is similar and omitted. □ 
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1.8.2 Proof of Theorem 4 



(Existence) For any h GTi and x, £ G M define the matrices ^f{x) by 



(1.115) 



The proof of existence rests on the martingale property of M+ = {M^}t>Q and M = {M^ }t>o given by 
M+ = e-fr^ ''(^^)<iV+ (^M./>^tA./) and Mf = e'/o"^' (^m.,- C*a.,- ) (l-HS) 



with 



/+(z,x) := e-^|(x)fc+, /-(t,a;) := e-^^ (a;)A;_, 



(1.117) 



where and fc_ arc N'^— and iV — column vectors, respectively. 
The martingale property of M+ follows from the equality 



where {Qt}t>a denotes the filtration generated by (^, Y). To verify this identity observe first that the Markov 
property of {^,Y) yields that 



(x 



Further, in view of the strong Markov property and spatial homogeneity of ^, the expectation on the right- 
hand side of the previous display is for a; < ^ given by 



= E, 



{Yo+=3,r^«^} 



E 



E ^0 

jeE°UE+ 

e-vt exp (Q+(f - x)) = e'^+{x)k+ = f+{i,x) 



e-x : 



(1.118) 



where the last line follows by the definitions (1.106) and (1.98) of and 77^. 

As M+ is a martingale, an application of Ito's lemma shows that /+ = {f+{i,u),i G E) satisfies for all 

u<(. 

\s{iff'l{i,u) + m{i)f+{i,u) + J2Qij{f+{j,u) - f+{i,u)) = 0, (1-119) 
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where and /" denote the first and second derivatives of /+ with respect to u. By substituting the 
expressions (1.115) - (1.117) into equation (1.119) we find, since fc"*" was arbitrary, that and r]^ satisfy 
the first set of equations of the system (1.66). The proof for and r]f^ is analogous and omitted. 

(Uniqueness) Now we turn to the proof of the uniqueness of the Wiener-Hopf factorization. To this end, 
let {W^ ,G~) be a complex matrix Wicncr-Hopf factorization and define the function / as in 
(1.117), but replacing r]+ and (5+ by and G+ respectively. Since the pair (1^+,G+) satisfies equation 
(1.66), it follows by an apphcation of Ito's lemma, that = e~ -fo ^^"^'^^^ f{Yt,^t) is a local martingale. In 
view of the facts that G+ € D(A''+) and h G H it follows that M' is in fact bounded on {t < r^}. An 
application of Doob's optional stopping theorem then yields that 



f{j,x) = E^j 



+ lim 



e-^«''<^^^'V(l^t,6)/{.+=J. (1.120) 



By the definition of /, the absence of positive jumps of ^ and (1.118), the first expectation in (1.120) is equal 
to 



:E., 



for X < i. The second term in (1.120) is zero, since Jq I{YseEa}^^ ^ oo P^j almost surely on the event 
{t^ = oo} for all j & E and x G M, and minjg^;,, 3?(/i(i)) > 0. Thus = f for all A^+-column vectors fc+ 
and we deduce that G+ = Q~l and = r/'^. Similarly, one can show that G~ = Qf^ and W~ = ri^ and 
the uniqueness is proved. □ 



1.8.3 Proof of Theorem 5 

Applying the strong Markov property at f and noting that f < yields, in view of the representation (1.94) 
that 



{Yf=j,T+<T^ } 



F{v,Yr) 



'J^HY,)ds 



{Yf=j,T+>T-} 



F{£,Yr) 



where 



since h{i) = h{i)I{i^EO}, by definition of h. From the definition of it is straightforward to check that 

= 0] = 1 ioriGE^ and x G {£, v}, 

so that F{v, i) = e"^5,j and F{i, i) = e"^Sij if i e E^, where Sij denotes the Kronecker delta. 

Moreover, in view of the form (1.61)-(1.62) of Qh and the definition of phase-type distribution, it is clear 
that, conditionally on Iq = * € Ej' and = Y^( = j and ^ PH(^j,B^), where Si is the vector with 
elements 6i = (Sik)- Therefore we find, using (1), that for u with Ji(u) < at and i G Ef 
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= e""E,,4e"^'] = e"" [{-ulf - Bf)-\-Bf)\\ {i). 
Similarly, it follows that for i G E~ and u with 3?(u) > —a~ 

F{t,i)=e^'[{uIJ -Bj)-\-Bj)l] ii). 
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